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The inert polymer polyethylene glycol (PEG) induces a "bundling" phenomenon in F-actin so- 
lutions when its concentration exceeds a critical onset value Co. Over a limited range of PEG 
molecular weight and ionic strength, C can be expressed as a function of these two variables. The 
process is reversible, but hysteresis is also observed in the dissolution of the bundles, with ionic 
strength having a large influence. Additional actin filaments are able to join previously formed 
bundles. PEG polymers are not incorporated into the actin bundles. 

Estimates of the Asakura-Oosawa depletion force, Coulomb repulsion, and van der Waals potential 
are combined in order to explain the bundling effect and hysteresis. Conjectures are presented 
concerning the apparent limit in bundle size. 

PACS numbers: 87.16.Ka, 82.35.Rs, 82.35.Pq, 61.25.Hq 



I. INTRODUCTION 

Cells and tissues are crowded with macromolecules 
such as proteins, DNA, and various polymers. While the 
ligand-substrate specific binding model forms the cen- 
tral dogma of biochemistry, many protein-protein inter- 
actions are facilitated in large part by physical drives 
such as the excluded volume effect, electrostatic inter- 
action, and depletion force. Chemically inert polymers 
such as polyethylene glycol (PEG) and dextran are often 
added into solutions of biomolecules in order to mimic the 
crowded biochemical environment and understand vari- 
ous biomolccular functions ranging from protein filament 
assembly , ion channel opening and closing , to tran- 
scription of DNA [3| ■ This study focuses on the physical 
effects of PEG on the lateral aggregation of a filamentous 
protein assembly, with the goal of a first principle expla- 
nation of its aggregation property. Similar phenomena 
occur in many cellular and physiological settings, which 
by large are dictated by the common physical mecha- 
nisms, although often in more complicated and less de- 
fined conditions. 

The protein selected for this study is actin. Actin 
is a ubiquitous cytoskeletal protein of molecular weight 
(MW) 42000 Dalton. In solutions of low ionic strength 
it exists as a globular monomer, G-actin. As [K + ] or 
[Na + ] is increased beyond 50 mM, the monomers poly- 
merize into helical filaments (37 nm pitch) of side-by-side 
monomers incorporated into two strands 0, 0, known 
as filamentous actin, or F-actin. F-actin is polydisperse 
in length, and consists of 370 monomers per /itm unit 
length. The diameter of F-actin is 8 nm. At pH 8, 



the linear charge density is about 4 e/nm and surface 
charge density 0.15 e/nm 2 . This value is derived from 
the amino acid sequence of a-skeletal muscle actin, where 
each monomer has a net charge of — 12e [4J, |6(. In a di- 
lute solution F-actin is a freely undulating filament with 
a bending modulus k c such that the persistence length 
L p = k c /kT = 17 iixsx 0, 13 This length is larger than 
an averaged filament length of F-actin. In the study of 
lateral aggregation of F-actin as induced by flexible poly- 
mers such as PEG, the effect of flexibility of F-actin is 
negligible, and the filaments are treated as charged col- 
loidal rods. 

The physical insight into the interaction of free poly- 
mers with colloidal particles was first achieved by 
Asakura and Oosawa (AO) a half century ago 9J. The 
essential prediction from the AO treatment is that an at- 
tractive force is generated between two colloidal particles 
in the presence of non-interacting polymers. The effect, 
known widely as the depletion force, is essentially of en- 
tropic origin, and has been calculated for various sizes 
and geometries of the colloids [HJ, El H E 13 

To apply the AO model, flexible polymers are typi- 
cally treated as freely interpenetrating hard spheres of 
radius Rao, which are excluded from the colloid sur- 
face by a thin layer of thickness Rao ■ It was shown 
that this shell creates a positive free energy difference, 
AF = PV = P Rao A c , where P is the osmotic pressure 
due to the polymer, and A c the surface area of the colloid. 
If two colloidal particles share part of this volume, the 
volume accessible to the polymers is increased. Conse- 
quently, the total entropy of the system is increased and 
thus the free energy of the system is lowered. In terms 
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FIG. 1: Cross-section of two aligned cylinders with a deple- 
tion^ exclusion, layer. Ra is the F-actin radius, R g is the 
PEG radius of gyration, and D is the axis-to-axis separation 
of the actin cylinders. On the left is a cartoon of the two 
cylinders surrounded by PEG, represented as self-avoiding 
random walks on a cubic lattice, N=180, corresponding to 
PEG 8000 [OJ. In this case C is well below C* . In the 
AO scheme on the right, the colloidal rods are surrounded by 
interpenetrating hard spheres of effective radius Rao- The 
attractive interaction is PSV, where 5V is the product of the 
hatched overlap area and an arbitrary length L. 



the subject of continuing investigation. Nevertheless, in 
a continuum approach Eqn.nbecomes a volume integral, 
AF = f P(r)dV. 

Two analytical treatments have been proposed recently 
to calculate the correlation length £ and osmotic pressure 
P in the semidilute regime. The renormalization group 
theory 0] yields expressions [l7L ITsj for £ and P as func- 
tions of polymer concentration and molecular weight. A 
second treatment has been developed by Schweizer and 
co-workers [Tol Il9j , which employs the polymer reference 
interaction site model (PRISM). The numerical results 
from these two treatments agree with each other in a 
range of concentrations spanning the dilute and semidi- 
lute regimes. The predictions also agree well with the 
available experimental data. Therefore, results from both 
treatments are applied later in this work in calculating 
the depletion force. 

Interacting charged colloids in saline solution are typ- 
ically modeled as bodies with Coulomb repulsion and 
attractive short-range dispersion forces, a combination 
known as the DLVO theory [53, [2l| . The electrostatic 
interaction in this situation is modeled with the Poisson- 
Boltzmann equation, which describes the relation be- 
tween charge density p and electrical potential <fi, under 
the assumption that the ionic charges deviate from their 
bulk concentration uq according to the Boltzmann law. 
For monovalent ions: 



of the Helmholtz free energy, 

AF(V,T) = ^6V = PSV. (1) 

For an ideal gas, ^ = ^, thus it is clear that the AO 
interaction is entropically driven. However, osmotic pres- 
sure and volume of exclusion are more accessible to mea- 
surement and calculation. The geometry of the colloids 
and their depletion layer is portrayed in Fig^ The vol- 
ume of the exclusion layer is a function of the polymer 
radius of gyration R g , the radius of the colloidal rods Ra, 
and the axis-to-axis separation between the rods D. It 
has been shown for an ideal random chain polymer 
that Rao = %Rgl V^, where R g is the radius of gyration. 



Polyethylene glycol (PEG) has long been used as an 
agent for precipitation and crystallization of proteins. 
PEG is a polymer that finds water a "good" solvent, 
in that it assumes the configuration of an almost ideal 
self-avoiding random walk |15| . with radius of gyration 
R g ~ MW 3 / 5 Hi'. As the polymer concentration C ex- 
ceeds a value C*, the molecules begin to overlap, and 
direct scaling arguments 0] show that the correlation 
length £ begins to decrease as (C/C*)~ 3 / 4 . The thick- 
ness of the polymer depletion layer at the colloid sur- 
face |l2j scales with £, with the polymer monomer num- 
ber density continuously passing from zero at the colloid 
surface to its bulk value. The nature of this profile is 



V 2 </> - -p/e 



p = e(n+ - n_) = n a e[e^ e ' kT - e^' kT ] 

or 

V 2 (f) = 2n o esinh(0e/A:T)/e. 

With boundary conditions of moderate surface charge, 
the solutions of this equation are associated with a char- 
acteristic exponential decay length known as the Debye 
screening length which for monovalent ions is pro- 
portional to 1/ ^/hq. 

Dispersion forces are the result of mutually induced 
dipole attraction between two bodies, and as such fol- 
low a 1/r 6 law for interaction potential [2(|, and are of 
course shorter in range than the electrostatic interaction. 
At high enough salt concentration, the electrostatic re- 
pulsion between colloidal particles of like charge can be- 
come weaker than the attractive dispersion force. As a 
result, the colloids precipitate. 

The rest of the paper is organized as follows. Section II 
describes materials and methods employed for this work. 
Section III presents quantitative data and microscopic 
observations about the formation of F-actin bundles that 
result when the PEG concentration exceeds a critical on- 
set value C . Much of the analytical effort is made in 
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section IV, by employing what are basically estimates 
of the electrostatic repulsion, van der Waals attractions 
(the DLVO theory), and the depletion effect. Qualita- 
tive explanations are made concerning onset of bundle 
formation, hysteresis, as well as the bundle size limit. 
Limitations of the present treatment and the experiments 
performed are also discussed. 



II. MATERIALS AND METHODS 

A. Materials 

G-actin was prepared from acetone powder of rabbit 
skeletal muscle after Spudich and Watt [22J and stored 
as 0.2 mL aliquots, 7.9 g/L, at -80 °C until use. G- 
actin was polymerized by first diluting to 3.2 g/L with 
"G-buffer", and then adding 3 M KC1 to bring [KC1] to 
150 mM. Polymerization to F-actin was immediately evi- 
dent by a rapid increase in viscosity. Polymerization was 
always allowed to proceed for several hours. "G-buffer" 
is 0.2 mM CaCl 2 , 0.5 mM ATP, 0.5 mM NaN 3 , 0.5 mM 
DTT, 2.0 mM tris-HCl at pH 8.0. To create stock solu- 
tions where [KC1] > 0.3 M, the salt was added directly to 
G-buffer. PEG of molecular weights 4000 to 20000 were 
purchased from Sigma and PEG 35000 from Alfa-Aesar. 
Radiolabeled PEG was obtained from American Radio- 
labeled Chemicals (ARC-1565, 3 H -labeled PEG 35000) 
and 14 C labeled PEG 4000 from Amersham Biosciences, 
Inc. (CFA508). 



B. Light Scattering Measurement of C a 

Ninety degree light scattering was measured by adopt- 
ing a Perkin-Elmer LS-5B luminescence spectrophotome- 
ter using a 1 mL culture tube (4.5 mm ID) to hold the 
sample. The tube was aligned with the probe beam for 
the least refraction by the glass tube and maximum illu- 
mination of the sample. Monochromatic light of 550 nm 
wavelength was used for illumination and detection, with 
spectral slit width selected from 3 to 15 nm to keep the 
detected light in the linear range of the photomultiplier. 
Once a measurement was commenced by adding the ini- 
tial PEG for a certain KC1 concentration, the spectropho- 
tometer settings were kept unchanged for consistency. 

Concentrated PEG solution was typically added in vol- 
umes of 5 /zL for a typical increment of 0.2% in the w/w 
percentage of PEG. Thorough mixing was achieved using 
a long, 0.2 mL pipet tip with a cut end (~ 1 mm dia), 
and shear was minimized by keeping the flow rate below 
0.2 mL/1 sec. During the first couple of mixing cycles, 
the contrast in the index of refraction in the mixture was 
seen to quickly diminish. This observation suggests that 
the mixture became homogeneous following several cy- 
cles of pipeting. F-actin concentration was 0.3 g/L with 
a sample volume of 0.4 mL. C Q was characterized by a 



sudden rise in turbidity, developing in less than 60 sec- 
onds, and a corresponding increase in light scattering. 



C. The amount of PEG in the bundle structure 

Approximately 2000 Bq (2000 counts/s) of 3 H labeled 
PEG 35000 was added to 0.7 mL F-buffer, 600 mM KC1, 
in a glass test tube of 1 cm diameter. The sample was 
then mixed with 0.2 mL of 3.2 g/L F-actin. To in- 
duce bundling, 10 /xL increments of 13% w/w unlabeled 
PEG 35000 in F-buffer containing 0.6 M KC1 were added 
with thorough yet gentle mixing. When the light scatter- 
ing (turbidity) characteristic of bundling was observed, 
one more 10 [iL aliquot was added, and the mixture was 
allowed to sit for 10 min. This resulting w/w concen- 
tration of PEG solution for the onset of bundling is de- 
noted C a . Then the entire sample was transfered to a 
1.5 mL plastic Eppendorf vial and centrifuged at 1200 g 
for 5 min. in a swinging bucket rotor. The supernate 
SO was extracted and the pellet was immediately resus- 
pended in 100 ^L of 7% w/w unlabeled PEG F-buffer 
containing 600 mM KC1. At this point 1.2 mL more of 
the same unlabeled PEG buffer was added and any re- 
maining clumps resuspended by pipet. To minimize shear 
in resuspension, the flow rate was limited to 100 /iL/sec 
in pipeting, and the pipet tip had the nozzle cut for a 
0.5 mm opening. 

After allowing the resuspension to re-equilibrate for 30 
min, the protein was pelleted again at 8000 g for 5 min. 
in a swinging bucket rotor, and the supernate SI ex- 
tracted. To minimize the pellet mass, the vial was again 
centrifuged at 13000 g for 5 min. in a fixed angle rotor, 
and the last remaining supernate removed with capillary 
action. The total weight of the vial was then obtained 
to 50 /ig accuracy, and the pellet weight W pc n t derived 
from the vial tare. The pellet was then suspended in 
100 fiL of 2% SDS and quantitatively washed into 5 mL 
of scintillation fluid (IRL BioSafe II). 10 /iL samples of 
SO and SI were added to 5 mL of scintillation fluid, and 
the activity measured with a Beckman LS6500 scintilla- 
tion counter. 



D. Hysteresis in bundle formation 

For a given KC1 concentration (100 mM or 600 mM) 
and 0.3 g/L F-actin, PEG 35000 was added in 0.2 % in- 
crements, with thorough but gentle mixing. When the 
qualitative change in turbidity associated with bundle 
formation was observed, the mixture was allowed to equi- 
librate for 30 minutes. Then, the sample was split into 
0.3 mL aliquots, and each centrifuged at 2000 g. An 
appropriate volume was removed from the top, and re- 
placed with F-buffer to reduce the PEG concentration. 
The bundles were then gently resuspended with a cut 
pipet tip. After 12 hours at 4 °C, the bundles were pel- 
leted at 2000 g x 7 min. and the supernate removed. 



4 



Centrifugation was repeated to remove remaining traces 
of F-buffcr. The bicinchoninic acid method (Sigma B- 
9643) was used to assay the remaining bundled protein 
in the pellet. DTT from the F-buffer was found to inter- 
fere with this assay, but was not significant in relation to 
the final amount of measured protein. 

E. Numerical solution of cylindrical PB equation 

A 'shooting' technique p3| , with potential <ft as the in- 
dependent variable [2J|, was applied to numerically gen- 
erate solutions to the PB equation with boundary con- 
ditions for an infinite cylinder. Exploiting the symme- 
try of the problem, we set 4>i(ri) = iA<\> for i = 0, 1, .... 
Debye-Huckel theory tells us that the counter charge den- 
sity p between and r^+i is 2n o esinh(0 i e/fcT)/e, where 
4> i = {4> i + <p i+1 )/2. Charge conservation is expressed 
with the relationship 

A 4 +i = A., + Ar2i:r i p((j) i ) 

where A^ represents the net linear charge density enclosed 
within the ith shell. This, with the following approxima- 
tion of Gauss's Law 

A<j> _ A i+ i 
Ar litTit 

allows a quadratic equation to be solved for Ar: 

A<f> = Ar(Ar2nr iP (fa) + X l )/2nr t e 

After an initial guess of ro = B, <fio = for the outer 
boundary condition, is calculated iteratively. If, at 
n ~ -Ractin, Ai — -A ac tin, the solution set n(4>i) is ac- 
cepted. If not, a new B is chosen. A value of 0.1 mV 
= kT/250 is used for A<f>. Typically the charge balance 
condition is satisfied to within 2 %. 



curves can be fit with the following semi-empirical func- 
tion: 

C o = a + ax/{l + [KCl]/5 mM) 2 (2) 

where C is percent w/w of PEG. a\ is found to be a 
rather strong function of molecular weight. A powerlaw 
fit yields ax = 70.0(MV4 7 /8000)" 2 - . Separate measure- 
ments at 600 mM KC1 show that clq is a weak function 
of molecular weight, a = I.0(MVy/8000) -0 ' 5 . The term 
1/(1 + [KC1J/5 mM) could be interpreted as the result 
of mass action ion association between the total bound 
charge of the F-actin Qq and the solution potassium ions 
of concentration [K + \. 

Ka 

Q~ + K+ ^ QK 



Hence 

[Q-] = [Q Q }/(l + K A [K + ]) 

where = 1/(5 mM) is the association constant in this 
ligand binding model j2(|. One can speculate that an 
effective charge [Q~] dominates the electrostatic interac- 
tion at lower ionic strength, where the repulsion between 
two F-actins, proportional to [Q~] 2 , is overcome by a de- 
pletion force proportional to the osmotic pressure at the 
onset of bundling. Similar empirical fits can be applied 
to the mobility of fd virus as a function of monovalent 
and divalent ionic strength [ Q. Wen and J.X. Tang, un- 
published data]. 



III. RESULTS 
A. Threshold and reversibility of bundle formation 

Four parameters of phase separation were explored: 
PEG MW, PEG concentration C, ionic strength, and 
actin concentration. The solution pH was held fixed at 
8.0. As polymer concentration C is increased, a sudden 
increase in light scattering marks the concentration C a 
where the onset of bundling is induced 0, |2^|. At 
monovalent ionic strength near 100 mM and PEG MW 
8000 or below, the turbidity of the sample is low as the 
structures are small. These small structures were also 
measured by fluorescent imaging (data not shown). For 
35 kD PEG, the light scattering is stronger, as the struc- 
tures are much larger in both length and diameter. 

The effects of PEG MW and the solution ionic strength 
on C are shown in Fig[3 The data clearly suggest an in- 
terplay of Debye screening length and polymer R g . These 



Measurements show small effect of protein concentra- 
tion on C over the range of 0.2 g/L to 0.8 g/L F-actin, 
although the kinetics for bundle formation, as indicated 
by light scattering, were appreciably slower for the lowest 
F-actin concentration (data not shown). 

Also, the bundling process is reversible. After PEG 
concentration passes C , the protein may be sedimented 
by centrifugation and resuspended in F-buffer with a re- 
sulting disappearance of turbidity. PEG may then again 
be added to induce bundles. 



Hysteresis is seen in the disappearance of the bundles 
if the polymer concentration is increased beyond C a and 
then, by dilution, decreased below C . As shown in Fig 
there is a pronounced effect of ionic strength, with 
bundles formed at a higher ionic strength more stable. 
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FIG. 2: The effects of PEG MW and ionic strength on 
bundling critical concentration C , as detected by the onset 
of light scattering. The lines through the points are semi- 
empirical fits using Eqn.|5| where Ka = 1/(5 mM) for each 
data set, but with an adjustable baseline oo for each MW. 
F-actii 
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FIG. 3: Hysteresis in the dissolution of PEG 35000-induced 
bundles at [KC1] = 100 mM and 600 mM. The PEG concen- 
tration is normalized to C/C , and the amount of recovered 
protein (bundles) is plotted as the fraction of the total pro- 
tein. For C < Co, there were virtually no bundles formed. 
As PEG is added, almost all the F-actin went into bundles 
and was pelleted. After this, as the PEG concentration was 
lowered in suspension, many of the bundles did not dissolve 
until C was lowered by over 50 %. 



induced by PEG 35000. It should be noted that there is 
no role of PEG concentration heterogeneity in the matu- 
ration, as the polymer concentration is entirely homoge- 
neous after the mixing process completed within the first 
few seconds. It was found, however, through separate ex- 
periments, that additional mixing at the later stages also 
affected the kinetics of the bundle growth. The physical 
mechanism of this additional observation is clear in that 
as the bundle size increases, diffusion becomes progres- 
sively inefficient in facilitating the self-assembly 

When F-actin is added to a solution where C = 5C D , 
large aggregates with sheet- or balloon-like structure are 
formed. This must largely be associated with the macro- 
scopic, slow mixing process of viscous media, easily seen 
as discontinuities of index of refraction on the mm length 
scale. Furthermore, actin will coalesce without mixing of 
the PEG solution. If a 10 /iL drop of F-actin is placed 
on the surface of a denser 5% PEG solution, the concen- 
trated protein is found as one mass on the surface. 

Under all circumstances, there were no signs of slow 
spinodal decomposition or coexistence of phases. 



It appears that individual labeled actin filaments show 
no size preference when they combine with previously 
formed, unlabeled bundles. When new phalloidin-labeled 
F-actin is introduced, the fluorescence seems to be uni- 
formly distributed. As is in Fig the DIC images co- 
localize with that of the TRITC fluorescence, suggest- 
ing that individual F-actin filaments attach themselves 
to the previously formed bundles. Also, it was possible 
to show that preformed bundles do not exchange mate- 
rial, as phalloidin-labeled bundles mixed with unlabeled 
bundles show separate DIC and fluorescence populations 
(data not shown). 



C. PEG is not part of the bundle structure 



B. Nucleation and growth 

Immediately after the polymer concentration reaches 
C , actin filaments coalesce into thin strings (FigEJ left). 
After a few minutes, with mild intermittent mixing, the 
typical size increases (Fig0J middle). Then, with mild 
intermittent mixing over 30 min, these initial structures 
grow into polydisperse structures larger in diameter and 
length (FigQ] right). The bundles in the first image have 
roughly 10 5 F-actin filaments each; those in the final im- 
age are at least an order of magnitude more massive. 
Therefore, some consideration must be given to the ki- 
netics of bundle formation. The series of images in Fig 
21 shows a typical time progression of the actin bundles 



As described in Methods, an elaborate procedure was 
designed to address the question whether PEG is an in- 
tegral part of the F-actin bundle structure. Following a 
systematic set of measurements, it was concluded that 
less than 0.1% w/w PEG was found to be tightly bound 
or trapped with the bundles. This result directly con- 
firms that PEG-induced bundle formation is the result 
of the depletion effect, rather than a putative PEG bind- 
ing or cross-linking of F-actin. 

IV. DISCUSSION 

We interpret the data in Fig Elas an interplay between 
the attractive depletion force and the electrostatic repul- 



FIG. 4: F-actin bundle formation and progression following the gradual addition of PEG 35000 to just past C . The three 
representative pictures show bundles immediately after the gradual addition of PEG 35000 to just past C (left panel), after 
5 min. (middle), and after 30 min. (right). [KC1] = 100 mM. Images are 160 x 160 fim, taken with a Nikon TE300 microscope, 
10X phase contrast. The samples were flowed into a channel of 70 /im thickness created by placing two pieces of double sided 
tape between a microscope slide and a cover slip. F-actin concentration is 0.2 g/L. 




FIG. 5: Left panel: 40X fluorescence image of PEG 35000 bundles mixed with sparsely labeled (TRITC-phalloidin) F-actin. 
Right panel: 40X DIC image at same location. The sample thickness in these images is less than 10 fj,m, so bundles are much 
more sparse in these images than in Fig [I] Image size is 170 /im x 135 /jm. 



sion of two negatively charged cylinders in a monovalent 
salt solution. We present approximations showing that 
it is reasonable to expect the depletion force to cause F- 
actin to laterally aggregate. Inclusion of a van der Waals 
attraction offers an explanation of the large hysteresis at 
high ionic strength, although it is not clear exactly how 
much this interaction contributes to the final bundled 
state. 

In the following three sections we calculate the electro- 
static, depletion, and van der Waals potentials between 
two infinite, parallel F-actin cylinders in order to un- 
derstand the onset of bundle formation. Based on the 
results of these calculations, issues such as the observed 



hysteresis and bundle size are discussed. 



A. Depletion Attraction 

As described in the Introduction, the depletion layer 
and osmotic pressure create an additional free energy for 
a single colloidal particle in a polymer solution. If the 
layer thickness were constant for all polymer concentra- 
tions and the osmotic pressure difference linear in concen- 
tration, this positive free energy would be J PdV, where 
P continuously passes from Pi, at the colloid particle sur- 
face to zero in the surrounding solution. Since this is not 
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the case, the energ y is found by integration over volume 
and concentration [L7l | 

JVJC' b =0 U ^b 

where C(r) is the local polymer concentration; the bulk 
polymer concentration is C' b , and the local osmotic pres- 
sure Pb is a function of C' b . This expression for integrating 
total free energy is in analogy with an electrostatic charg- 
ing process to calculate, for example, the energy required 
to charge a capacitor. 

Following the qualitative reasoning of the AO theory, 
the free energy for two interacting colloidal particles, as 
expressed in Eqn. should become less positive as the 
two approach, indicating an attractive depletion force. 
To calculate the interaction, one places the two particles 
at separation D and zero polymer concentration. Then 
the polymer bulk concentration is parametrically ramped 
to the final value Cb, with Eqn. 0] used to evaluate W. 
This process has the qualitative and intuitively appealing 
feature of generating a depletion force with a range of the 
order of R g for all concentrations. We denote this change 
in free energy W as the particles approach from infinity 
as Ud- 

We employ osmotic pressure data and £ values from 
the PRISM theory |10l , and an ansatz about overlapping 
depletion regions [Ljj to estimate the attractive depletion 
potential Ud- For osmotic pressure P(C), we adopt a 
virial equation P = B\C + B2C 2 from the literature |10| . 
The values for B\ and B2 thus determined are presented 
in Table HI 

We use the proposal 01 that the polymer concentra- 
tion in the overlap region C(r) can be expressed as a 
product of the separate concentrations: 

d(r) C a (r) 
C(r) = — j^Cb (5) 

where Cb is the bulk concentration and C\ and C2 are 
the profiles surrounding the two colloids at infinite sep- 
aration. C(r) is in effect the product of two probabil- 
ity densities. The concentration profile near an isolated 
cylinder is assigned a cone-shaped function of width 2 £ : 

r — R a 

Ci{r) = <3>— — A R A <r<R A + 2i 

With the virial expressions for P(Cb) and the deple- 
tion profile C(r) from Eqn. [3] we numerically evaluate 
Eqn. Letting Cb take the values of C in Fig [21 for 
PEG 35000, namely 1.6, 1.08, 0.97, and 0.80 % w/w (cor- 
responding to 100 mM, 140 mM, 175 mM, and 500 mM 
KC1 respectively), we show the numerical results of Ud 
in Fig This more elaborate methods of calculating 
the attractive potential Ud yield values less than the AO 
hard sphere model, typically by a factor of five. 



TABLE I: A list of 4 parameters for PEG used in this study. 
B\ and B2 are the first and second viral coefficients in the 
formula P = B\C + B2C 2 where C is percent w/w and the 
osmotic pressure P is in the unit of Pa. The values for R g and 
C* are derived from Kulkanari et. al. flTjl and are consistent 
with other measurements [l5|. 



PEG MW 


Bi 


B 2 




C*(%w/w) 


8000 


3112.5 


771.1 


4.7 


3.0 


20000 


1245.0 


641 


8.2 


1.4 


35000 


711 


574 


11.4 


0.9 



B. van der Waals attraction 

The standard formula for the van der Waals interac- 
tion |2l| for two parallel cylinders, based on the Der- 
jaguin approximation, seemed to have an unreasonably 
long range, so a numerical integration of the potential 
was done, based on the 1/r 7 force between two elements 
of identical material. Assuming two elements of material 
separated by r experience a force F = A-^d 3 rid 3 r2/rl 2 
where A is the Hamaker constant [2£j , we see that an el- 
ement dx\dy\dz\ and a rod of area dx2dy2 separated by 
S experience a force F = A-^j^dxidyidzxdxzdyz/ S 6 . 
Considering a circular slice of F-actin, height dz, and a 
parallel cylinder of F-actin with lattice size dz, one can 
numerically calculate the total F for a series of sepa- 
rations D, then numerically integrate Fdx to obtain a 
potential U v dw(D) for this slice. This function is shown 
in Fig ©for A = 0.1 kT, L = 100 ran, 

C. Electrostatic Repulsion 

For a single cylinder of surface charge density a, in 
salt solution of ionic strength /, the potential solution of 
the linearized Poisson-Boltzmann equation is known to 

be mm 

<j}{r) = aKo{Kr)leKK 1 {R A K) (6) 

where Kg and K\ are modified Bessel functions, and k _1 
is the Debye screening length [2(|; for monovalent salts, 
kT x = 0.3 nm/vl where I is the value of ionic stength 
in moles/L. The numerical results described in Methods 
match closely Eqn. [HI for I = 100 mM. The inset of Fig 
[H] shows numerical solutions <fi(r) for KC1 from 100 to 
600 mM. 

We should mention in passing that with divalent coun- 
terions, and for cylinders of smaller diameter and higher 
surface charge, (f>(r) does not have the approximate form 
4> — </>oexp(— (r — Ra)k). Under these conditions there 
is a sudden rise in potential as r — > Ra- For a cylin- 
der of modest surface charge, the system free energy is 
4>oQ/2 m,m|' I n this situation, the surface potential 
is less than kT/e = 25mV, and the positive entropic 
energy of the counterion cloud is equal to its negative 
electrostatic energy 30] , so change in system free energy 
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FIG. 6: Numerical values for U v dw, Ud, and Ue- For Ue, 
monovalent salt ranges from 100 mM to 600 mM. For the Ud 
family, PEG 35000 concentration has the values in Fig 
1.60, 1.08, 0.95, and 0.80 percent w/w. The cylinder radius 
and surface charge are those of fully charged F-actin, R — 
4 nm, A = 4 e/nm, For U v dw, a Hamaker constant A — 0.1 kT 
was chosen. L = lOOnm. Inset: (j>{r) for monovalent salt from 
100 mM to 600 mM. 



is reflected as change in surface potential. For F-actin in 
100 mM KC1, the electrostatic potential at the surface is 
calculated to be 30 mV (see inset of FigEJ. Therefore, F- 
actin with monovalent counter ions is in the linear regime 
of the PB equation. In contrast, DNA (5.9 e/nm, 0.9 nm 
diameter) has a significantly higher surface potential, and 
non-linear effects must be taken into account. 

We do not have a calculation of electrostatic interac- 
tion Ue{D) for two parallel cylinders in proximity. How- 
ever, because the highest value for (f>(r) (the surface po- 
tential) is not much greater than kT/e, we use the sum 
(pi (r) + <f>2 (t) as an approximate solution: 



V 2 (</>!+ 02 ) 



. , ,<P\e . foe 

sinhl ) 

y kTkT' 



smh(— )- 



-smh(— ) 



We take the change in free energy as the two cylinders 
are forced together to be like two charges in free space: 



U E (D) 



i(T\dS\ 



(7) 



Si 



For D > 2R + this is a very reasonable approxi- 
mation [2^. We neglect the effect in the region where 
the dielectric of one cylinder displaces the counterions 
of the other, assuming the protein cylinder has the 
same dielectric constant as the surrounding medium. In 
Fig [5] we show the numerical integration of Eqn. {7\ for 
[KC1] = 100 mM to 600 mM for fully charged F-actin. 
These results match the low surface charge example of 
Harries El- 



10kT 



-lOkT 



- ■ lOOmM KC1, 1.6% p35kD 

— 500mM KC1, 0.8% p35kD 
■■■ 500mMKCl,0% p35kD 




10 



15 



D (nm) 



20 



25 



FIG. 7: Sum of Ue, Ud, and f7 v dw for fully charged (4e/nm) 
F-actin. Values for [KC1] and PEG-35kD are taken from the 
data of FigH L = 100 nm. 



D. Prediction of bundle formation 

Our interpretation is that the depletion potential is 
sufficient to allow the filaments to find a lower energy 
condensed phase, which is normally denied by Coulomb 
repulsion, a view consistent with the standard DLVO the- 
ory. The qualitative feature to be appreciated is the role 
of ionic strength in making the vdW interaction accessi- 
ble. 



Shown in Fig is the sum Ud + Ue + U v aw for 
three representative conditions: 500 mM KC1, 0.8% 
PEG 35kD; 100 mM KC1, 1.6% PEG 35kD; and 500 mM 
KC1, 0% PEG. The calculations appear successful in ex- 
plaining features of the results at 500 mM ionic strength, 
predicting bundle formation with 0.8% PEG 35kD, con- 
sistent with Fig [21 With a rather low Hamaker constant, 
A = 0.1 kT, the electrostatic repulsion is large enough 
to stabilize the suspension at 0% PEG, as it must. Ac- 
cording to Eqn. [7| fully charged F-actin in 100 mM KC1 
has a large value of Ue that totally prohibits the role of 
van der Waals forces in the bound state. This is contra- 
dictory with the experimental result. Both 100 mM and 
500 mM PEG 35kD bundles are similar in appearance, 
although hysteresis has a wider range for 500 mM KC1, as 
seen in Fig [21 While the choice of the low Hamaker con- 
stant A = 0.1 kT is an obvious cause of the imbalance, it 
is hard to justify using different values for Hamaker con- 
stant at different ionic strengths. Additional mechanisms 
unaccounted for in our treatment are discussed below. 

Counterion fluctuation and redistribution [2U l3ll I32I 
l33l | clearly play important roles in leading to an over- 
all attractive interaction. These effects have not been 
taken into consideration in our simple electrostatic cal- 
culations. Since the effects are more pronounced at lower 
ionic strength, it is not surprising that our calculated Ue 
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based on PB theory overestimates the repulsion (Fig |BJ). 
at low ionic strength. Most of these works seem to focus 
on polyelectrolytes of high surface charge density such as 
DNA, in the presence of multivalent counterions. These 
theories may also be expanded to treat systems of rela- 
tively weak surface charge, such as F-actin. 

The semi-quantitative fits to the data shown in Fig |3 
suggest ligand binding as an alternative mechanism for 
predicting bundle formation. With a mass action lig- 
and (counterion) binding interpretation (Eqn. |3J) with 
Ka= 1/(5 mM), half the negative charge would be neu- 
tralized at 5 mM KC1, which implies a surface charge 
down by a factor of 100 at 500 mM. This surface charge 
could certainly not support a stable suspension. A 
weaker association constant, say Ka = 1/(100 mM), fails 
to produce the flat baseline of C at higher salt. Ex- 
periments |3 focused mainly on F-actin's tightly 
bound divalent ions Ca 2+ and Mg 2+ have demonstrated 
effects of [K + ] in the rapid phase of fluorescence induc- 
tion, indicating low affinity K + binding constants from 
1/(10 mM) to 1/(100 mM). The data in Fig together 
with the fact that F-actin is stable at the highest mono- 
valent ionic strength, seem to demand that F-actin main- 
tain a significant fraction of unneutralized charge at high 
salt. Indeed, the electrophoretic mobility measurements 
for filamentous phages fd and M13 suggest that the fila- 
ments remain charged up to 500 mM KC1 [Q. Wen and 
J.X. Tang, unpublished data]. 

A recent study by Yu and Carlsson j^| directly ad- 
dresses the question of F-actin electrostatic interaction. 
This work includes electrostatic and entropic terms in 
ligand (counterion) binding as well as the actual spacial 
configuration of the charged groups of the protein struc- 
ture. These authors point out a role of induced charge 
condensation as the two filaments begin to interact: as 
the surface potential of one filament is increased by the 
other, the effective Ka of each site is raised by a fac- 
tor exp(e02 / kT) , lowering the repulsion by increasing 
the counterion binding. The actual interaction is found 
with a largest-error-correction algorithm seeking the low- 
est energy, which is a function of the charge of all sites. 
Their treatment may be expanded to predict features of 
PEG-induccd bundle formation, which goes beyond the 
scope of this work. 



E. Distinction from polyvalent cation induced 
bundle formation 

The mechanism of polymer-induced bundling is dis- 
tinct from the bundling induced by divalent metal ions, 
cationic complexes, and basic polypeptides |2^, l3^ |. 
The most salient feature is the inhibiting nature of mono- 
valent salts for the latter vs. the enhancing effect for the 
former. Secondly, polyvalent counterions such as basic 
polypeptides and cationic complexes function as some 
binding agent in F-actin bundles. An opposite role is 
held by PEG, which does not bind F-actin, as confirmed 



in this study with the sedimentation assay using radio- 
labeled PEG. In addition, we show in this work that un- 
der some conditions the PEG-induced bundles are mor- 
phologically different in size. 

Another distinct property from polyvalent cation in- 
duced bundle formation is that no rcsolubilization was 
found with high polymer concentration or high ionic 
strength (2.5 M KC1), in that the structures were sta- 
ble under these conditions. In contrast, for most col- 
loidal systems, precipitation occurs with the addition 
of multivalent ions, only to go back into suspension 
with increased concentration of the multivalent counte- 
rions. Long known for some classical colloids, this has 
recently been shown for the bacteriophages fd and M13 
viruses |37j . This phenomenon has been predicted for 
F-actin by recent theoretical calculations |36j. However, 
experimental test of such a prediction has yet to be per- 
formed. 

There is an apparent multitude of proteins which ei- 
ther bind to F-actin, affect the polymerization of actin, 
or enhance the lateral aggregation of F-actin [25j. In 
view of the distinction between the polyvalent counte- 
rion induced bundling and the one by depletion effect 
shown in this work, it is helpful to assess their respec- 
tive contributions to the related biochemical functions. 
Under the unifying tie of electrostatic effects, however, 
this work implies a strong relevance of polyelectrolyte 
physics to the observation that so many biochemically 
unrelated proteins bind F-actin and induce formation of 
actin bundles. Is the binding capability due largely to 
the cylindrical geometry of actin filaments, which tends 
to diminish electrostatic stability [2?ll30| ? 

F. Hysteresis 

The hysteresis in the dissolution of PEG induced actin 
bundles (Fig 3) can be qualitatively explained by the 
following scenario: At zero polymer, filament/filament 
contact is prohibited by the electrostatic barrier. As the 
height of this barrier is overcome with the depletion po- 
tential Ud from increased polymer concentration, F-actin 
is allowed into the vdW binding configuration. If the 
polymer concentration is then lowered, the tightly bound 
state remains stranded inside the Ue barrier. This can 
be seen in Fig [7| for the case of 500 mM KC1. 



G. What limits the bundle size? 

Suppose that the total surface charge of a bundle in- 
creases with size, while the surface charge density re- 
mains constant. Recall that Eqn. shows that the self- 
energy goes as X 2 /VlR ~ a 2 R/VT. If the surface charge 
density is a constant value (that of an individual F-actin), 
this electrostatic self-energy per unit length of a bundle 
would increase as the bundle radius R. Assuming the 
depletion potential to be of depth £, Ud per unit length 
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would scale as £ x / 2 ii 3 / 2 , which is the shaded area of Fig 
n for small £/R. /.From these considerations, there is 
no obvious limit in size because Ud outraces Ue with 
increasing diameter 2R. This simple scaling argument 
also predicts that linear polyelectrolytes of larger diam- 
eter are more prone to bundling when other parameters 
are comparable. 

The final bundle structure is not necessarily a hexag- 
onal array of rods, but rather it could be the product 
of diffusion limited aggregation (3^ . Such structures are 
expected to show surface roughness greater than the stan- 
dard statistical deviation proportional to the square root 
of the accumulation. From Fig [5] we know that an in- 
dividual filament finds a low energy configuration on a 
previously formed bundle. However, the situation is dif- 
ferent for two matured bundles, where, due to surface 
contour, the relative interaction (contact) area fails to 
scale with diameter. Such a structure is easily broken up 
by the turbulence of light mixing. Thus there is a point 
in the hierarchy of assembly when the macroscopic hy- 
drodynamic forces exceed the surface-surface interaction. 
This argument explains why the bundles do not collapse 
into one large mass, but assume polydisperse structures 
which show no tendency to combine with each other after 
a certain level of self-assembly. 

The arguments above based on the surface roughness 
do not fully explain how bundles might reach some equi- 
librium size limit. The experimental observations are in- 
conclusive whether an equilibrium size distribution exists 
or is practically attainable. PEG 8kD bundles appear to 
stop growing at smaller sizes (data not shown) than those 
induced by PEG 35kD (Fig |j| and Fig , which clearly 
grow into large structures, their extent perhaps only lim- 
ited by the progressively slow kinetics. In the test tube 



experiments, the effect of diffusion may be surpassed by 
the sedimentation effect due to gravity for the large bun- 
dles observed. Additional experiments are necessary in 
order to assess which effects are dominant in determining 
the actin bundle size. 

The hysteresis data also implies an inhomogeneity, or 
variation in 'fitness', among the actin bundles. Some of 
the structures 'survive' with reduced polymer while oth- 
ers have dissolved back to dispersed actin filaments. The 
inhomogeneity may be related to defects in the struc- 
ture of the bundles, or more simply a wide variation of 
the bundle sizes. Additional experiments may also be 
designed to test these different proposals. 



H. Concluding remarks 

Our calculations are to be taken as reasonable esti- 
mates of the forces involved in bundle formation. More 
exact calculations, especially of the electrostatic interac- 
tion, presumably including a role of mass action counter- 
charge accumulation, should help explain many aspects 
of the results quantitatively. Nevertheless, results from 
the simple experiments reported here provide an oppor- 
tunity to expand the understanding of some basic prob- 
lems in lateral aggregation of protein filaments. 
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